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ABSTRACT 



Context. The magnetic and kinematic properties of the photospheric Evershed flow are relatively well known, but we are still far from 
a complete understanding of its nature. The evolution of the flow with time, which is mainly due to appearance of velocity packets 
called Evershed clouds (ECs), may provide information to further constrain its origin. 

Aims. We undertake a detailed analysis of the evolution of the Evershed flow by studying the properties of ECs. In this first paper we 
determine the sizes, proper motions, location in the penumbra, and frequency of appearance of ECs, as well as their typical Doppler 
velocities, linear and circular polarization signals. Stokes V area asymmetries, and continuum intensities. 

Methods. High-cadence, high-resolution, full vector spectropolarimetric measurements in visible and infrared lines are used to char- 
acterize the ECs phenomenon through a simple line-parameter analysis as well as more complex Stokes inversions. 
Results. We find that ECs appear in the mid penumbra and propage outward along filaments with large linear polarization signals and 
enhanced Evershed flows. The frequency of appearance of ECs varies between 15 and 40 minutes in diff'erent filaments. ECs exhibit 
the largest Doppler velocities and linear-to-circular polarization ratios of the whole penumbra. In addition, lines formed deeper in the 
atmosphere show larger Doppler velocities, much in the same way as the "quiescent" Evershed flow. According to our observations, 
ECs can be classified in two groups: type I ECs, which vanish in the outer penumbra, and type II ECs, which cross the outer penumbral 
boundary and enter the sunspot moat. Most of the observed ECs belong to type I. On average, type II ECs can be detected as velocity 
structures outside of the spot for only about 14 min. Their proper motions in the moat are significantly reduced with respect to the 
ones they had in the penumbra. 
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1. Introduction 

The Evershed flow has been s tudied intensiv ely over the last 
' decades (for a review see, eg. ISolankil[20Q3h . However, most 
of what we know today about this dynamical phenomenon has 
been learned from observations that do not reflect its temporal 
I variability. The study of the evolution of the flow requires sta- 
ble observing conditions during long periods of time, which is 
difl&cult without the help of adaptive optics systems. For this 
reason, ther e are very f ew investigation s published about the 
topic (Shine et al.lll994l: Ikimmele 1 99j;^ Rouppe van der VoortI 
l2QQ3l: [Georgakil as & Christopouloull2QQ3l) . Their main result is 
that the Evershed flow displays variations on time scales of 8-25 
min. This quasi-periodic behavior is produced by the so-called 
Evershed clouds (ECs), velocity packets which propagate from 
the mid penumbra to the outer penumbral boundary. 

A good knowledge of ECs is crucial for a complete under- 
standing of the nature of the Evershed flow. Unfortunately, their 
observational properties are not well characterized yet. Previous 
investigations were based on high resolution filtergrams and 
Dopplergrams, but the polarization signals of ECs have never 
been observed. As a result, important details such as the relation 
between magnetic fields and ECs remain unknown. 
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In a series of two papers we present the first study on the 
temporal evolution of the Evershed flow that makes use of high- 
cadence, high-precision, full vector spectropolarimetric observa- 
tions of visible and near-infrared lines. The spatial resolution of 
these measurements, around Off 6 at 1565 nm and Offl at 630 nm, 
is one of the highest ever achieved in solar spectropolarimetry. 

This first paper concentrates on the observational charac- 
terization of ECs. In Sect. [2] we describe the observations and 
data reduction. We then identify the ECs and determine their 
lifetimes, proper motions, and sizes (Sect. [3]). Sections |4] and 
\5\ investigate the properties of the penumbral filaments along 
which ECs move, and the periodicity of the EC phenomenon. 
We find that ECs occur in the center of intra- spines, i.e., elon- 
gated structures having the mor e horizontal field s and stronger 
Evershed flows of the penumbra (iLites et al.ll 19931) . The spectro- 
scopic and polarimetric properties of the ECs are presented in 
Sect. [6l and compared with those of the intra-spines in Sect. [71 
The dependence of the Doppler velocities on height is examined 
in Sect.m We also study the correlation between proper motions 
and Doppler velocities (Sect. O, the evolution of the properties 
of the ECs as they propagate outward (Sect.[TO|), and their disap- 
pearance inside (Sect. [TT]) and outside (Sect. [TJ]) the penumbra. 
Our main findings are summarized in Sect. \13\ 

Paper II of the series ('Cabrera S olana et al.l2QQ7 l) contains an 
analysis of the thermal, magnetic, and dynamic properties of the 
ECs as inferred from inversions of the measured Stok es spectra. 
Further details can be found in . Cabrera Solanal (l2QQ7l) . 
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Table 1. Visible and infrared spectral lines observed with TIP 
and POLIS. Aq represents the laboratory central wavelength,;^ 
the excitation potential of the lower level, and log gf the loga- 
rithm of the oscillator st rength times the n i ultiplicity of the lowe r 
leve l, a, b, and c refer tofNaye et alJ (ll994 . lBorrero et al.l(l2003l) . 
and iPierce & Breckenridgd (Il974l) . respectively. Note that O2 
are telluric absorption lines. 
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Fig. 1. Speckle-reconstructed G-band filtergram of AR 10781 
taken at the DOT on 2005 June 30, 9:34 UT. The box shows the 
fraction of the center- side penumbra and adjacent moat scanned 
by TIP and POLIS. The arrow marks the direction to disk center. 



Table 2. Log of the VTT observations. The second column gives 
the heliocentric angle of the spot. The third and fourth columns 
show the scanning step and the integration time per slit position. 
The duration of the sequences and the seeing conditions are in- 
dicated in the last two columns. 



Date 


6 


step size 


h 


duration 


seeing 




[deg] 


[arcsec] 


[s] 


[min] 


[arcsec] 


06/30/2005 


43 


0.2 


10.5 


156 


~ 0.6-0.7 


07/01/2005 


35 


0.2 


5.5 


80 


~ 0.8-0.9 



2. Observations and data analysis 

2.1. Observations 

NOAA Active Region 10781 was observed on 2005 June 30 
from 8:47 to 1 1 :23 UT and on 2005 July 1 from 9:3 1 to 10:5 1 UT 
at the German Vacuum Tower Telescop e (VTT) on Tenerife. The 
Tenerife Infrared Po larimeter (TIP: iMartmez Pillet et all 1 19991: 
ICollados et al.lll999 | ) and the Polarimetric Littrow Spectrograph 
(POLIS: ISchmidt et al.l[2QQ3l: iBeck et al.l[2005bh were operated 
simultaneously to record the full Stokes profiles of the spectral 
lines around 1565 nm and 630 nm (Table [T]). The same spot 
was observed on 2005 June 30 with the Dutch Open Telescope 
(DOT) from 8:45 to 9:38 UT. Table |2] gives details of the obser- 
vations at the VTT. The active region was located at heliocentric 
angles of 43° and 35° on June 30 and July 1, respectively. The 
slit width was a.'36 for TIP and ^^8 for POLIS. During the 
observations, the Kiepenheuer Adaptive Optics System (KAOS; 
ISoltau et aDl2002 ) was used to reduce image motion and blur- 
ring. In order to observe the same field of view (FOV) the TIP 
and POLIS slits were aligned by rotating the main spectrograph 
of the VTT and moving the TIP camera. The pixel size was 
(TfllS and a.' 145 for TIP and POLIS, respectively. 

On June 30, we performed rapid scans of a small portion 
of the center- side penumbra of the spot including the adjacent 
moat region (see Fig. [T]). The scan step was Qf!2 for a total of 
20 steps (4''). The integration time was 10.5 s per slit position, 
resulting in a cadence of 3.9 min. We performed 40 repetitions of 
the scan. The excellent seeing conditions, together with KAOS, 
allowed us to reach a spatial resolution of about Of! 6 at 1565 nm 
and Qffl at 630 nm. These values have been derived from the 
power spectra of the quiet sun granulation surrounding the spot. 



On July 1 we observed a smaller portion of the center- side 
penumbra of the same active region: the scan step was Of! 2 for 
a total of 10 steps, producing maps of T'. The integration time 
was 5.5 s per slit position (2 min cadence). The seeing condi- 
tions were good and the spatial resolution was around QffS in the 
infrared and Off 9 in the visible. 

The use of infrared (IR) and visible (VIS) spectral lines, with 
diff'erent heights of formation, allows us to obtain physical infor- 
mation from a wide range of layers in the atmosphere. Moreover, 
the diff'erent sensitivities of the spectral lines to the atmospheric 
parameters impr ove the diagnostic capabilit ies of a single spec- 
tral range alone (ICabrera Solana et al ] l2QQ5h . 

2.2. Data reduction 

2.2.1 . Calibration of the polarimetric data 

To remove the diff'erent response of the CCD pixels to homo- 
geneous illumination, gain tables calculated using flat field and 
dark current images were applied. Then, the intensity signals 
measured by the polarimeters can be converted to Stokes profiles 
using the demodulation matrices derived with the help of the in- 
strumental calibration unit of the VTT. Instrument al crosstalk 
was re moved applying the telescope model of iBeck et aP 
(l2QQ5al) . while residual crosstalk due to the coelostat mirrors was 
corrected for using the statistical method of ICollados! (l2QQ3h and 
Schlichenmaier & Collados ( 2002) . The accuracy of the polar- 
ization profiles is estima ted to be on the ord er of 10"^ in units of 
the continuum intensity (iBeck et al .112005 ah . 

2.2.2. Scattered light correction 

The presence of residual O2 telluric lines in the polarization 
profiles recorded by POLIS reveals the existence of spectrally 
undispersed scattered light (Fig. IA.2I) . If a constant fraction of 
the incident light enters the camera without going through the 
grating, the intensity measurements are contaminated by a cer- 
tain amount of light artificially produced by the instrument. This 
scattered light might be due to incomplete shielding of the light 
path inside the compact POLIS shelter. 

We correct the spectra f or scattered light in the same way as 
lAllende Prieto et aL l (l2Q04l) . The amount of scattered light and 
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Table 3. Mean dispersions of the Doppler shifts of the selected 
umbral pixels for the different spectral lines. The errors represent 
the standard deviations of the averages. 



Date 


O' 630.15 


O' 630.25 


O' 630.35 


O' 1564.85 


O' 1565.29 




[m s-i] 


[m s-1] 


[m s-1] 


[m s-i] 


[m s-i] 


06/30 


83 + 16 


80 + 17 


100 + 20 


86 + 16 


120 + 16 


07/01 


64 + 30 


63 + 31 


74 + 34 


64 + 28 


60 + 28 



the resolving power of the polarimeters are calculated by com- 
paring their average quiet-sun intensity profiles with those in the 
Fourier Transform Spectrometer (FTS) atlas. The FTS profiles 
are not aff'ected by scattered light and the point spread function 
(PSF) of the FTS can be considered as a Dirac delta function 
for our purposes. Hence, if we add a certain fraction of scattered 
light to the FTS profiles and convolve them with the spectral PSF 
of the polarimeters, we should be able to reproduce the observed 
line shapes, i.e., 

where / are the observed disk-center profiles, /fts the FTS pro- 
files, K is the fraction of scattered light, G the polarimeter in- 
strumental profile (assumed to be a Gaussian of width cr), and ^ 
represents convolution. 

Figure lA.ll displays an example of the fits obtained us- 
ing Eq. [H The best results are achieved with K ^ 15% and 
cr ^ 6 pm, both for the July 1 and June 30 data sets. The the- 
oretical resolving power of POLIS (2 pm, limited by the slit) is 
always smaller than that inferred from the fits. The somewhat 
large value of cr may be a consequence of imperfect focussing 
of the instrument. 

2.3. Wavelength scale 

The velocity scale for the IR an d VIS lines has been set assuming 
that the umbra does not move (lBeckerslll977|) . The wavelength 
scale is defined by the rest wavelength and the spectral disper- 
sion. We obtain the former calculating the mean zero-crossing 
positions of Stokes V for the pixels located in the darkest regions 
of the umbrcQ. For the June 30 time sequence, we select pixels 
with continuum intensities ///qs < 0.4 at 630 nm. The disper- 
sion is computed by determining the number of pixels between 
the line core positions of Fei 1564.85 nm and Fei 1565.52 nm 
for TIP and the O2 telluric lines for POLIS. 

To estimate a lower limit for the calibration error, we cal- 
culate the dispersion of the Doppler shifts in the selected um- 
bral pixels with respect to their mean value for each map. In 
Table [3] we show the average of those quantities in the two se- 
quences. The errors indicate the standard deviations of the aver- 
ages. As can be seen, the mean dispersions are always smaller 
than 120 m s"^ with low scatter (rms < 3 4 m s~^) betw een maps. 
Similar dispersions have been found by iRezaei et"aD (12006) in 
their analysis of velocities in sunspot umbrae. Such small val- 
ues could be produced by th e non-zero velocities of umbral dots 
(ISocas-Navarro et aLll2004l) and/or the presence of oscillations 
in the selected areas. 

^ For the POLIS data, the small focal length of the spectrograph pro- 
duce a curvature of the spectrum along the slit. In order to remove it, 
we use the O2 telluric line at 630.20 nm to align all the spectra individ- 
ually. This correction has been applied before the determination of the 
wavelength scale. 
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Fig. 2. From left to right: Maps of speckle-reconstructed G-band 
intensity, continuum intensity at 1565 nm, and continuum inten- 
sity at 630 nm observed on June 30, 9:26 UT. Contours outline 
pixels with G-band intensities greater than 0.5 /qg. Each tickmark 
on the }7-axis represents V\ 

2.4. Alignment of TIP, POLIS, anci DOT data 

2.4.1 . Alignment of TIP and POLIS 

Since the refractive index of the air depends on wavelength, the 
images of the Sun taken by T IP and POLIS are shifted with re- 
spect to each other (R eardonll2QQ6|) . In o rder to remove th ese 
shifts, we follow the method described by lBeck et al.l (l2007h . It 
basically consists in performing a cross-correlation of the inten- 
sity maps using the largest pixel (the TIP pixel) as a reference. 
To that aim, the POLIS data are rebinned to the TIP pixel size. 
Figure [2] displays, as an example, continuum intensity maps at 
1565 nm and 630 nm after removing the shifts. As can be seen, 
the correspondence between the intensity structures is quite pre- 
cise. The ac curacy of the alig nment between TIP and POLIS is 
around Of! 1 (iBeck et al.ll2QQ7l) . 

The position of the Sun on the sky changes during the day. 
Hence, the diff'erential refraction produces time-dependent shifts 
between the VIS and IR spectropolarimetric data. When the time 
sequences are long, one has to take into account the temporal de- 
pendence of the displacements. This is why we calculate the dis- 
placements between IR and VIS spectra for each repetition of the 
scan. The total displacements accumulated during the 156 min of 
observations on June 30 are around V and 0'.'6 in the directions 
perpendicular and parallel to the slit, respectively. This means 
that CO- spatial Stokes profiles in the visible and the near-infrared 
are not strictly simultaneous, but the time diff'erence is always 
smaller than ~ 55 s. 

2.4.2. Alignment of VTT and DOT data 

The way filtergrams and spectropolarimetric maps are obtained 
is diff'erent. The DOT filtergrams are 2D images taken at once, 
while the maps recorded at the VTT are constructed by scan- 
ning the slit across the solar surface, which takes time. If one 
wants to compare both kind of observations, the positions ob- 
served by the polarimeters have to be locate d on the DOT fil - 
tergrams taken at the same time. Following iBeck et al.l (l2007l) . 
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Fig. 3. Evolution of line parameters in the center-side penumbra of AR 10781 on June 30. Top to bottom: Continuum intensity at 
1565 nm, Stokes V zero-crossing velocity of Eei 1565.28 nm and Eei 630.25 nm, logarithm of the total linear polarization of Eei 
1564.85 nm and Eei 630.25 nm. Color contours outline the ECs. The letters at the bottom of the third panel label each EC. Red 
lines indicate the inner and outer penumbral boundaries. Each tickmark in the j-axis represent V\ The arrow marks the direction to 
disk center, t = min corresponds to June 30, 8:47 UT. 



we align the POLIS slit-jaw images with the DOT G-band im- 
ages. This is done by degrading the G-band data to the pixel size 
of the slit-jaw images. After the rotation between the images of 
the two telescopes is removed, shifts between the slit-jaw and 
the G-band images are determined by cross-correlation. At this 
point, one knows where the positions of the POLIS/TIP slits are 
in the G-band images. Einally, one constructs artificial G-band 
maps by selecting the co-spatial G-band strip which is closer in 
time to the polarimetric data. The DOT and VTT observations 
are not strictly co-temporal, but the differences are small enough 
(< 30 s) to work with both data sets as if they had been taken 
simultaneously. 



2.5. Line parameters 

The simplest way to analyze the observations is by computing 
line parameters from the polarization profiles. Those parameters 
reflect the atmospheric conditions of the solar photosphere, pro- 



viding a first estimate of their values. We have derived the fol- 
lowing observables from the Stokes profiles of the IR and VIS 
spectral lines: 

1. Continuum intensities, /. According to the Eddintong- 
Barbier approximation, / depends on the temperature of 
the layers where the continuum is formed [~ log(T) = 
01. This approximation works well in sunspot pen umbrae 
(Idel Toro Iniesta et al.lll994l:lBellot Rubio et al.|[2Q06) . 

2. Stokes V zero-crossing velocity. The Doppler velocity is es- 
timated from the Doppler shift of the Stokes- V zero-crossing 
wavelength, Azc- ^zc is calculated as the mid-point between 
the maxima of the Stokes V lobes. This velocity is represen- 
tative of the magnetic atmosphere inside the resolution ele- 
ment and can only be computed for pixels exhibiting normal 
two-lobed Stokes V profiles. 

3. Total polarization signals. Eor fully split lines, if the mag- 
netic field or ientation remains constant along the LOS 
(ISolankilll993h : 
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Fig. 4. Same as Fig.E from t = 78.4 to t = 156.6 min. 



- The total linear polarization, L = Yji ^Q{M)^ + U(Ai)'^, 
is proportional to sin^yLOS- Hence, the larger the value 
of L, the greater the inclination of the magnetic field to 
the LOS (rLos). 

- The total circular polarization, V = 2/ is propor- 
tional to C0S7L0S- Therefore it is smaller when the field 
inclination is larger. 

- The linear-to-circular polarization ratio is greater when 
the inclination increases (L/V ~ sin^ Tlos/ cosyLOs)- 

The polarization signals provide only a rough estimate of 
Tlos. because the magnetic field is known to vary with height 
in the solar photosphere and the spectral lines are not always 
in the strong field regime. 
4. The Stokes V area asymmetry, given by 



SA 



\Ar\ - lAbl 
|Ar| + |Ab|' 



(2) 



where and A^ are the areas of the red and blue 
lobes of Stokes V. SA arises as a consequence of ve- 
locity gradients along the LO S (lAuer & HeaslevI Il978l: 
ISanchez Almeida & LitesI [T992I) . The combination of such 



gradients with magnetic field or temperature gradients leads 
to an amplification of the asymmetries. 
We construct magnetograms, M, as (minus) the Stokes V sig- 
nal of Fei 630.25 nm at AA = -\-lO pm from line center. M 
provides a rough estimate o f the longitudinal magnetic flux 
(iLandi degrinnocentill992h . 



3. Identification, morphology, and proper motions 
of ECs 

In Figs. [3] and m we show the temporal evolution of the center- 
side penumbra and the moat region of AR 10781 during the 156 
min of the June 30 data set. The first panel displays the con- 
tinuum intensity at 1565 nm. The next two panels display the 
Doppler velocities calculated from the Stokes V zero-crossing 
shifts of Fe 1 1565.28 nm and Fe 1 630.25 nrr0. The last two pan- 
els show the logarithm of the total linear polarization of Fei 
1564.85 nm and Fei 630.25 nm. Appendix |B] gives the same 
information for the July 1 data set. 



^ The velocity is computed only for pixels exhibiting normal two- 
lobed V profiles with amplitudes greater than 0.2% for Fei 630.25 nm 
and 0.5% for Fei 1564.85 nm. 
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Table 4. Phenomenological properties of the observed ECs. ECs reaching the quiet photosphere are marked in bold face. Parentheses 
indicate ECs that leave or appear outside of the FOV. The date and time of appearance of each EC are specified in the second and 
third columns. The fourth column displays the radial distance at which the EC is seen for the first time. Only the values for the ECs 
appearing inside the FOV are given. The fifth column displays the time each EC stays inside the penumbra. The seventh column 
shows the propagation velocity, corrected for the viewing angle. The errors are the standard deviations in the determination of the 
proper motions. The last two columns display the maximum lengths and widths attained by each EC, corrected for projections 
eff'ects. 
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We identify ECs as structures of enhanced velocity signal 
which propagate outward across the penumbra. During the 236 
minutes of our observations, fifteen ECs appear in the mid- 
penumbra and then migrate to the outer penumbral boundary. 
They are observed in the two spectral ranges simultaneously. We 
have outlined them with contours of diff'erent colors and letters 
in Figs. [3] and m The boundaries of the ECs have been defined 
as isocontours of specific values of the Doppler velocity, hence 
they should be considered as approximate. 

In order to quantify the radial distances where the ECs ap- 
pear, we calculate the temporal average of the penumbral borders 
and define the mean penumbral radial distance, r, as the normal- 
ized distance from the inner (r = 0) to the outer (r = 1) mean 
edges of the penumbra. 

Table |4] summarizes some phenomenological properties of 
the observed ECs. The radial distances where they appear, rapp, 
vary from 0.4 to 0.7 with an average of 0.6. Their lifetimes, r, 
inside the penumbra range from 3.9 to 39 min with an average 
of 13.5 mir0. ECs can be classified in two diff'erent groups: (a) 
those that disappear at the outer penumbral boundary, hereafter 
type I ECs; and (b) those that reach the quiet photosphere and 
enter the moat (D, F, M, and N), type II in the following. Most 
of the ECs belong to the first group. 

The propagation velocities of the ECs inside the penum- 
bra, corrected for the viewing angle, range from 1.4 km s"^ 
to 4.2 km s"^ with an average of 2.6 km s"^ (seventh col- 
umn of Table |4l). These v elocities are in rough agreement 
with values reported earlier (Shine et al.lll994l: lRimmelelll994l: 
iGeorgakilas & Christopoulou 12003 ). To compute the propaga- 
tion velocity of each EC we have fitted straight lines to the (x, 3;) 
positions of the edge farther from the umbra as a function of 
time. The slopes of these lines give the apparent propagation ve- 
locities along the line of symmetr}0 (LS) and perpendicularly to 



it. Due to projection eff'ects the distances are contracted along 
the LS. We have corrected for this assuming that ECs move par- 
allel to the solar surface. Note that the vertical motion of ECs 
must be very small, otherwise they would leave the line forming 
region in a few minute^. 

In our observations, most ECs appear as elongated structures 
along the radial direction, and their sizes usually increase as they 
migrate outwards. To quantify the typical sizes of the ECs we de- 
fine a major and a minor axis for each one. The major axis is the 
segment connecting the edges of the EC along its longer direc- 
tion. The minor axis is perpendicular to the major axis and passes 
through its midpoint. We compute the length and width of the 
ECs along the two axes and correct them for projection eff'ects. 
The maximum length of the ECs varies from 500 to 3000 km 
with an average of 1700 km, while the maximum width ranges 
from 400 to 1200 km with an average of 700 km (cf. the two last 
columns of Table |4]). Note that for these computations we only 
use the ECs which are not partially outside the FOV. In Fig. [5] we 
show histograms of the maximum, minimum, and mean sizes 
along the major and minor axes for each EC. This figure also 
de monstrates th at the lengths are usually greater than the widths. 

IShine et al.l (119941) described the ECs as cloud-like struc- 
tures with size roughly 1000 km in length and a radial spac- 
ing of about 2000 km betwe en subsequent ECs . Our lengths are 
slightly larger than those of IShine et al.l (11994 "), while the typi- 
cal separation between two consecutive ECs is around 500 km 
in our observations (cf. second panels of Fig. [3] at t = 15.6 and 
t = 54.8 min). 



4. Relation with penumbral filaments 

The second, third, and fourth panels of Figs. [3] and |4] demon- 
strate that ECs move outward following penumbral filaments 



^ T is only calculated for the ECs that do not leave the FOV 

^ The line of symmetry is the projection of the LOS onto the solar 

surface. In the case of sunspots, it coincides with the line connecting 

the sunspot center and the disk center. 



^ The thickness of the photosphere is considered to be around 
500 km. If the vertical motions were greater than 600 m/s, ECs depart- 
ing from continuum forming layers would cross the entire photosphere 
in less than their mean lifetime (~ 13.5 min). 
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Fig. 6. From left to right. Maps of Stokes V zero-crossing velocity for Fei 1565.28 nm (top) and Fei 630.25 nm (bottom), total 
linear polarization, total circular polarization, and L/V ratio for Fe 1 1564.85 nm (top) and Fe i 630.25 nm (bottom) at ^ = 93.8 min. 
Color contours mark, and numbers label, three different penumbral filaments. White contours outline ECs F and G. The green line 
represents a cut along filament #2. 



with larger Doppler velocities and total linear polarization sig- 
nals than their surroundings. At least three diff'erent penumbral 
filaments can be distinguished in the June 30 data set. One of 
them is visible in the middle of the maps during the whole time 
sequence. The other two are partially outside of the FOV. We la- 
bel the left-side, center, and right-side filaments as #1, #2, and 
#3, respectively. The sixth column of Table |4] identifies the fila- 
ment along which each EC propagates outward. 

Figure [6] blows up the velocity and polarization maps of the 
penumbra at ^ = 93.8 min. The various panels display Doppler 
velocities, total linear polarizations, and total circular polariza- 
tions, along with their ratio, L/V. Penumbral filaments are out- 
lined with contours of diff'erent colors which have been con- 
structed taking into account the velocity and linear polarization 
signal maps. As can be seen, the filaments show up as struc- 
tures with enhanced v, L, and L/V, and their contrast is much 
higher in the quantities derived from the infrared lines. They ex- 
hibit larger circular polarization signals than their surroundings 
in Fe i 1564.85 nm, but smaller signals in Fe i 630.25 nm. These 
properties indicate that the filaments ob served in the maps are 
penumbral intra- spines (iLites et al.ll 19931) . 

In Fig. [7] the total linear polarization of Fei 1564.85 nm is 
shown as a function of time for the maps where EC A is de- 
tected. Pink and yellow contours outline EC A and filament #2, 
respectively. While the shape of the filament changes with time, 
it always remains visible from the inner to the outer penumbral 
boundary. We note that the pixels with larger values of the lin- 
ear polarization (marked with orange contours) coincide with the 



positions of EC A. This behavior is also observed in the other 
ECs. Hence, in addition to higher Doppler velocities, ECs dis- 
play larger linear polarization signals than the intra- spines along 
which they move. 

The polarization signals provide a rough estimate of the 
magnetic field inclination under some restrictive assumptions: 
if the pixel is occupied by a laterally homogeneous atmosphere, 
the lines are fully split, and the orientation of the magnetic field 
remains constant along the LOS, a higher inclination to the LOS 
implies larger L and smaller V signals, i.e., larger L/V ratios. 
Thus, Fig. [6] suggest that the intra- spines have more inclined 
magnetic fields than their surroundings, while Fig. [7] could in- 
dicate that the ECs are the most inclined structures inside them. 
However, the diff'erent behavior of the observables (in particular 
Stokes V) emphasizes the need of inversions for an unambiguous 
interpretation of the measurements. 



5. Periodicity 

As already mentioned, the Evershed flow displays variations on 
time scales of 8-25 min. The quasi-periodic behavior is pro- 
duced by ECs moving radially outwards along intra-spines. The 
perio dic times scales o f 15 min are associated with t he appear- 
ance dShine et al.ll 19941) or appearance/disappearance (iRimmeld 
199?) of ECs in isolated filaments. A precise estimation of the 
frequency of this repetitive behavior requires a Fourier anal- 
ysis of the data sets. Looking at Figs. [3]-[4l it is obvious that 
the orientation of the filaments changes with time. Hence, one 
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Table 5. Maximum Doppler shifts, linear-to-circular polarization ratios, and Stokes V asymmetries attained by each EC. The second 
and third columns are the unsigned blueshifts of Ee 1 1564.85 nm and Ee 1 630.25 nm. The fourth and fifth columns display the ratios 
of total linear and circular polarization of Eei 1564.85 nm and Eei 630.25 nm. The last two columns are the area asymmetries of 
Eei 1564.85 nm and Eei 630.25 nm. The errors represent the standard deviation of the individual parameters within each EC. The 
mean values of the observational parameters for the June 30 and July 1 data sets are marked in bold face. 



EC 


|Vl564.85l 


1^630.25 1 


(^/V)i564.85 


(^/V)630.25 


1564.85 


^^^630.25 




[km s-1] 


[km s-1] 






[%] 


[%] 


A 


3.3 + 0.3 


2.4 + 0.2 


0.44 + 0.02 


0.19 + 0.03 


5.9 + 2.0 


2.4 + 2.3 


B 


3.4 + 0.3 


2.5 + 0.3 


0.50 + 0.05 


0.27 + 0.02 


7.3 + 3.2 


7.7 + 1.8 


c 


3.2 + 0.3 


2.6 + 0.3 


0.45 + 0.03 


0.20 + 0.03 


6.9 + 1.7 


4.2 + 3.3 


D 


3.2 + 0.5 


2.4 + 0.5 


0.44 + 0.07 


0.22 + 0.04 


5.8 + 2.4 


5.1+2.8 


E 


3.1+0.3 


2.1+0.4 


0.44 + 0.06 


0.27 + 0.02 


5.3 + 1.9 


7.0 + 3.7 


F 


3.4 + 0.5 


2.5 + 0.4 


0.56 + 0.14 


0.29 + 0.05 


6.5 + 2.4 


9.6 + 6.5 


G 


3.0 + 0.3 


1.5 + 0.1 


0.55 + 0.05 


0.27 + 0.03 


9.1+3.7 


10.1 + 1.6 


H 


3.1 + 0.2 


1.9 + 0.3 


0.49 + 0.02 


0.24 + 0.01 


6.1 + 1.7 


5.8 + 2.5 


I 


3.1 + 0.2 


1.9 + 0.2 


0.36 + 0.04 


0.14 + 0.02 


3.6 + 1.3 


-2.4 + 1.3 


J 


3.2 + 0.3 


2.3 + 0.3 


0.49 + 0.02 


0.24 + 0.01 


5.5 + 2.6 


1.7 + 2.5 


K 


3.1+0.2 


1.5 + 0.4 


0.41 + 0.05 


0.16 + 0.02 


3.2 + 1.5 


-3.2 + 2.3 


Mean 


3.2 


2.2 


0.47 


0.22 


5.9 


4.4 


L 


2.1+0.1 


1.4 + 0.1 


0.87 + 0.07 


0.37 + 0.03 


4.7 + 1.0 


-2.1 + 1.1 


M 


2.2 + 0.1 


1.5 + 0.1 


0.83 + 0.06 


0.41 + 0.04 


5.5 + 1.4 


4.7 + 3.4 


N 


2.1+0.2 


1.2 + 0.2 


0.83 + 0.08 


0.38 + 0.03 


5.9 + 2.5 


6.9 + 4.7 


O 


1.9 + 0.3 


1.0 + 0.1 


0.69 + 0.02 


0.31+0.02 


5.4 + 2.7 


3.3 + 1.8 


Mean 


2.1 


1.3 


0.81 


0.37 


5.4 


3.2 




BCDEFGHI J 
Evershed cloud 




ABCDEFGHI JKLMNO 
Evershed cloud 

Fig. 5. Maximum, mean, and minimum lengths (top) and widths 
(bottom) of the ECs, corrected for projection eflfects (light to 
dark gray). 

should carefully select diff'erent paths for each map, otherwise 
the power spectrum will be contaminated with information from 
spatial points outside the filament and will show artificial peaks. 
Given its difl&culties, we do not pursue this kind of analysis here. 
However, a visual inspection of the maps provides valuable in- 
formation about the periodic appearance of ECs. 

In Eigs.[3]and|4]one can observe a repetitive behavior associ- 
ated with the appearance of ECs along filament #2. Its time scale 
is about 15 min. A cycle begins with EC A at ^ = 3.9 min. This 




0.0 



Fig. 7. Logarithm of the total linear polarization of Eei 
1564.85 nm as a function of time. Yellow and pink contours 
mark filament #2 and EC A, respectively. Orange contours out- 
line pixels with logLi564.85 > 1. 



EC arrives at the penumbral border around 15 minutes later. At 
that moment, another EC appears in the mid-penumbra (EC B) 
and a new cycle begins. After EC B reaches the outer penumbral 
edge, EC C appears and another cycle starts. New cycles take 
place in the filament with ECs D and E. This behavior is irreg- 
ular and does not happen during the whole time sequence. Eor 
filament #1 the periodic time scale seems to be of the order of 
40 min. This is the time gap between the appearance of ECs E, 
G, and I (third column of Table |4l). Eilament #3 is visible dur- 
ing the second half of the sequence. Only two ECs are observed 
to follow this filament, and the time between their appearances 
is around 23 min. Therefore, our observations suggest the exis- 
tence of diff'erent time scales associated with diff'erent filaments. 
Longer time series mapping larger regions of the penumbra are 
required to confirm or discard this possibility. 

6. Spectropolarimetric properties of ECs 

In this Section we determine the typical Doppler velocities, L/V 
ratios, and Stokes V slycsl asymmetries associated with ECs. To 
simplify the analysis, a single value of the parameters is calcu- 
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lated for each EC as follows. We compute the spatial average 
of the line parameters in each of the maps where a given EC is 
observed, and then take the maximum of these average values. 

Table \5\ shows the results for the spectral lines with higher 
sensitivity to magnetic fields (Fei 1564.85 nm and Fei 630.25 
nm). The maximum unsigned blueshifts are about 3.2 km s"^ and 
2.2 km s"^ for Fei 1564.85 nm and Fei 630.25 nm in the June 
30 time sequence, and 2.1 km s"^ and 1.3 km s"^ in the 1 July 
data set. These ve lociti es are lo\yer than the values reported by 
IShine et al.1 (1 19941) and lRimmell (Il994l) . The reason is that we 
calculate an average over a group of pixels, while they used ve- 
locities from individual pixels. 

The mean linear-to-circular polarization ratios are of the or- 
der of 0.47 and 0.22 in Fei 1564.85 nm and Fei 630.25 nm for 
the June 30 observations, and around 0.81 and 0.37 for the July 1 
data set. Since ECs show much larger L/V ratios in the infrared, 
vector magnetographs working at infrared wavelengths would 
detect them more easily than visible magnetographs. 

The Stokes V area asymmetry of Fei 1564.85 nm, M1564.85, 
amounts to 5.9% on June 30 and 5.4% on July 1, while M630.25 is 
around 4.4% and 3.2% on June 30 and July 1, respectively. For 
some ECs, the Fei 630.25 nm area asymmetry is of the order 
of its standard deviation, suggesting that ^530.25 has a strong 
dependence on radial distance. 

The large variations of the Doppler shifts and L/V ratios of 
the ECs as observed on two consecutive days are consistent with 
the idea that the flow and magnetic field are almost horizontal 
to the solar surface. The sunspot was nearer to disk center on 
July 1 . If the velocity and magnetic field vectors are horizontal 
to the solar surface, their projections to the LOS would cause 
smaller blueshifts, larger linear polarization signals, and smaller 
circular polarization signals (i.e., larger L/V ratios) as the spot 
approaches the disk center. This explains why the velocities are 
higher on June 30 while the L/V ratios are larger on July 1. 

7. Comparison of line parameters in ECs and 
intra-spines 

In a sunspot penumbra, the physical properties are known to vary 
radially from the umbra to the outer penumbral boundary. For 
this reason, we compare the observational properties of the ECs 
and those of the intra-spines as follows. We calculate the mean 
line parameters at each radial distance for the pixels that form the 
ECs and divide them by the mean parameters of the pixels inside 
the intra-spines which are at the same radial distance. This cal- 
culation is done for each map separately. Figures [8] and [9] show 
histograms with the results of the analysis for the Fe 1 1564.8 nm 
and Fe 1 630.25 nm lines observed on June 30. The results for the 
July 1 data set are equivalent (not shown). 

Interestingly, ECs exhibit larger Doppler velocities and L/V 
ratios than the penumbral intra-spines in both Fei 630.25 nm 
and Fe 1 1564.85 nm. Hence, our measurements suggest that ECs 
possess the most inclined magnetic fields of the penumbra. In 
Paper II we will confirm or disprove this conjecture by inverting 
the observed spectra. 

ECs also show larger Fe 1 1564.85 nm area asymmetries than 
the intra-spines. This might be the result of larger gradients of 
the physical quantities along the LOS and/or enhanced visibility 
of some magnetic components of the penumbra in the ECs. By 
contrast, M630.25 is smaller than the intra-spine average inside 
the ECs. This opposite behavior is surprising and could be due 
to the diff'erent heights of formation of the two spectral lines. 

The relation between ECs and continuum i ntensi ty is 
not well established in the literature. IShine et aP ([1991 " and 



Li 564.85/'^ 1564.85 




Fig. 8. Normalized distribution of the ratio of Doppler velocities 
(v), linear-to-circular polarization signals (L/V), Stokes V area 
asymmetries (SA), and continuum intensities (/) of the ECs and 
the rest of pixels of the intra-spines located at the same radial 
distance. The line parameters are derived from Fe 1 1564.85 nm. 
The vertical solid lines represent the mean of the distributions. 
The histograms correspond to the June 30 data set. 
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Fig. 9. Same as Fig.[8]for Fei 630.25 nm. 



iRouppe van der Voori (l2003l) report that ECs are mai nly cor- 
related with locally bright continuum features, while Rimmelg 
(1994) and Georgakilas & Christopoulou ( 2003) do not find any 
correlation. Just by visual inspection of the observations it is dif- 
ficult to discriminate if ECs are bright or dark continuum fea- 
tures. In the intensity histogram of Fig.O the ECs seem to show 
a preference for locally dark continuum intensities at 630 nm. A 
similar trend is observed in the continuum at 1565 nm (Fig. [8]), 
but less pronounced. These results have to be taken with caution 
because of the existence of unresolved structures at the angu- 
lar resolution of our observations (~ 0.6'' - 0.7''). In Fig. \TU\ 
we display co-spatial maps of speckle-reconstructed G-band in- 
tensity, continuum intensity at 630 nm, and continuum inten- 
sity at 1565 nm. The figure demonstrates that what we detect as 
dark/bright features in the spectropolarimetric data are indeed a 
mixture of both in the G-band filtergrams. As the ECs migrate 
outward, the dark/bright features inside them change, but in gen- 
eral the G-band intensities are locally darker at the position of 
the ECs. Spectropolarimetric observations at higher spatial reso- 
lution are therefore required to draw definite conclusions on the 
brightness of ECs. 

Finally, we note that the histograms of Figs. [8] and [9] pro vide 
information about the average line parameters of the ECs inde- 
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Fig. 10. From top to bottom: Cospatial speckle-reconstructed G- 
band filtergrams, continuum intensity maps at 630 nm, and con- 
tinuum intensity maps at 1565 nm. Color contours outline the 
ECs. 



pendently of their position in the penumbra. These properti es can 
chang e with radial distance. For example, Schlichenmaie r et al.l 
(l2QQ5l) suggested that flow filaments link bright and dark contin- 
uum features at diff'erent radial distances. Since ECs move along 
intra- spines (i.e., flow filaments), we may expect them to be as- 
sociated with locally bright features in the inner penumbra. This 
conjecture will be examined in Sect. [TOl 



3.5 
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Fig. 11. Maximum Doppler velocities of the observed ECs as 
derived from the Stokes V zero-crossing wavelengths. Maximum 
blueshifts for Fei 1564.85 nm, Fei 1565.28 nm, Tii 630.38 nm, 
and Fe i 630.25 nm are represented with diff'erent shades of gray, 
from light to dark. Note that the blueshifts are unsigned. 



9. Proper motions vs Doppler velocities 

It is of interest to compare the proper motions of the ECs along 
the LOS (calculated as the motion along the LS multiplied by 
the sine of the heliocentric angle) with their maximum Doppler 
velocities. If the blueshifts associated with ECs in the center-side 
penumbra were produced by proper motions, then both quanti- 
ties should be similar in magnitude. 

Figure [121 plots the LOS projection of the propagation ve- 
locities of the ECs vs the maximum Doppler velocities inferred 
from Fei 1564.85 nm and Fei 630.25 nm. The horizontal bars 
indicate the standard deviation of the Doppler velocities within 
each contour, and the vertical bars the errors in the determination 
of propagation velocities. The observed LOS velocities are sys- 
tematically higher than the LOS componen t of the proper mo- 
tion. This result agrees with the findings of IShine et all (flool 
and Rimmele (1994) and suggests that the Evershed flow can- 
not only be the manifestation of the radial propagation of ECs 
across the penumbra. The observed behavior i s also compati- 
ble wi th the sea-serpent scenario described by 'S chlichenmaied 
ioo2), in which magnetic flux loops slowly migrate outward 
with high flow velocity inside them. 



8. Dependence of Doppler velocities on height 

The magnitude of the Doppler velocity depends on the spec- 
tral line chosen to calculate it (cf. the first two columns of 
Table [5]). Figure [TT] shows the maximum unsigned Doppler 
velocities of the ECs inferred from each spectral line. There 
we see that the velocities indicated by the infrared lines 
are the largest ones. The IR lines of the set are formed 
deeper than the others. This implies that the Evershed flow 
is more conspicuous in deep atmospheric layers, in agree- 
ment with previous findings (iMaltbvl 11964: Balt hasar et al.l 
1 9971: IW es tendorp Plaza et al. l '200 ll: iMathew et aP 120031: 
Schfichenmaier et al.1 120041: Bellot R ubioetal.1 l2006h . 



The 



diff'erent velocities inferred from the diff'erent spectral lines may 
be the result of (a) a change of the modulus of the velocity 
vector with height; (b) a change of the inclination of the velocity 
vector with height; or (c) a combination of both eff'ects. 



10. Evolution of line parameters with radial distance 

In Sect. [6] we compared the line parameters of ECs and intra- 
spines disregarding the fact that there is a radial variation of the 
physical conditions in the penumbra. Here we study how the line 
parameters change with radial distance. To that aim, we average 
the parameters of pixels inside the ECs which are at the same ra- 
dial distance, and construct a curve of their variation with r. We 
note that those curves represent the average evolution of the ECs 
in the penumbra. To compare the behavior of the ECs and intra- 
spines, we do the same calculation for pixels inside the intra- 
spines excluding those of the ECs. 

Figures [13] and [141 show the results of this analysis for Fe i 
1564.85 nm and Fei 630.25 nm, respectively. As can be seen, 
the Doppler velocities are greater inside the ECs at all radial 
distances. Also interesting is the diff'erent shape of the velocity 
curves. The velocity always increases within the ECs as they mi- 
grate outward, from V630.25 = -1-8 km s"^ and V1564.85 = -2.8 
km s"^ at r = 0.5 to V630.25 = -2.2 km s"^ and V1564.85 = -3.0 
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Fig. 12. Proper motions along the LOS vs maximum Doppler 
velocities from Fe 1 1565.28 nm (left) and Fe i 630.25 nm (right) 
for the different ECs. Circles represent the maximum average 
Doppler velocities. The horizontal bars represent their standard 
deviations. The vertical bars are the errors in the determination 
of proper motions. 



km s"^ at r = 1. In the intra-spines, however, the Doppler veloc- 
ity slightly decreases from r ~ 0.6 outward. 

The ratio of linear-to-circular polarization is larger inside the 
ECs at all distances, which is more obvious for the IR line. In 
addition, both the ECs and the intra-spines display larger L/V 
values at larger radial distances. For the ECs, the L/V ratio in- 
creases from (L/y)630 25 = 0.14 at r = 0.5 and (L/y)i564 85 = 
0.38 up to (L/y)630.25 = 0.24 and (L/y)i564.85 = 0.47 at r = 1. 

The area asymmetry SA of the two lines also increases mono- 
tonically with radial distance, showing larger values within the 
ECs than in the intra-spines. The area asymmetries change from 
^^630.25 = -1.5% and Mi564.85 = 4.5% at r = 0.5 to M630.25 = 
6.2% and M1564.85 = 7.2% at r = 1. 

The continuum intensity curves display an interesting varia- 
tion with radial distance. Both /jr and /yis are greater than the 
intra-spine average inside the ECs at the radial distances where 
they appear (r ~ 0.4). From the mid to the outer penumbra, how- 
ever, the continuum intensities of ECs and intra-spines do not 
exhibit significantly diff'erent values. 

1 1 . Disappearance of type I ECs 

Most of the type I ECs observed on June 30 disappear suddenly 
in the penumbra, near the outer sunspot boundary. This indicates 
that the time they take to vanish is usually shorter than the ca- 
dence of the observations (3.9 min on June 30). There is, how- 
ever, an interesting case in which a type I EC is seen to disappear. 

Figure [15] shows the evolution of the Doppler velocities of 
EC A derived from Fei 1565.28 nm. As can be seen, EC A has 
vanished at ^ = 23.4 min (there is no dark patch where one would 
expect to find it, just outside of the penumbra). In order to quan- 
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Fig. 13. Average Stokes V zero-crossing velocity, ratio of the to- 
tal linear and circular polarization, area asymmetry, and contin- 
uum intensity of Fe i 1564.85 nm (top) for pixels inside the ECs 
(solid line) and the rest of pixels of the intra-spines (dashed line) 
as a function of radial distance, for the time sequence observed 
on June 30. Shaded areas indicate the standard deviations around 
the mean value. 



tify the evolution of the velocities in the EC, we choose radial 
cuts going from the inner to the outer penumbral boundary pass- 
ing through the middle of the EC in each map. These cuts are 
indicated by the vertical solid lines in Fig. \T5\ The variation of 
the Doppler velocities along the cuts is shown in Fig. [161 Colors 
indicate diff'erent times and the horizontal lines mark the upper 
and lower edges of the EC. In agreement with the behavior indi- 
cated by Figs. [13] and [141 the EC velocities increase with radial 
distance until t = 15.6 min. Then, there is a strong drop in the 
velocities at ^ = 19.5 min, after which the EC is no longer seen. 
This behavior is also observed for EC C and must be general for 
type I ECs. 



12. Evolution of type II ECs beyond the outer 
penumbral boundary 



IShine et al ] (119941) andiRi mmeld (1 19941) reported that some ECs 
penetrate slightly beyond the outer penumbral boundary before 
dissolving. On June 30 and July 1, four ECs crossed the edge of 
AR 10781. Here we describe the evolution of these type II ECs 
in the sunspot moat. 
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Fig. 16. Variation of |v63o.25l and |vi565.28l along the radial cuts shown in Fig.[T5l The colors indicate different times and the horizontal 
lines mark the upper and lower edges of EC A. 




Fig. 17. Temporal evolution of line parameters along intra-spine #2 (cf. the green radial cut in the first panel of Eig. O. Top: Stokes 
V zero-crossing velocity of Eei 1565.28 nm (left), continuum intensity at 1565 nm (center), and total linear polarization of Eei 
1564.85 nm (right). Bottom: Stokes V zero-crossing velocity of Eei 630.25 nm (left), continuum intensity at 630.2 nm (center), and 
total linear polarization of Eei 630.25 nm (right). The x-axis gives the average penumbral radial distance and the y-dod^ displays 
time. The arrows indicate the positions of ECs D and E. 



12.1. Modification of the outer edge of the penumbra by ECs 



In Eigure[6]we presented line parameter maps for filament #2. 
ECs D and E moved along this intra-spine. Eigure[T7]displays the 
temporal evolution of the Doppler velocity, the continuum inten- 
sity, and the total linear polarization as a function of the mean ra- 
dial distance for a cut along the intra-spine. The three signals are 
strongly correlated near the outer penumbral boundary, which 
becomes more evident when an EC reaches the quiet photo- 
sphere. This happens at ^ = 58.7 min for EC D and t = 91.1 min 
for EC E. When an EC penetrates into the quiet photosphere, the 
outer penumbra does the same: the velocity and linear polariza- 
tion increase and the continuum intensity decreases in the quiet 
photosphere, i.e., the penumbra gro ws at that position Such a 
behavior has already been reported by lShine et al.l (Il994|) for the 
intensity and velocity signals. 



12.2. Disappearance of type II ECs In the moat 

The evolution of ECs D and E outside of the spot can be followed 
in Eigs. 5 and 6. EC D (pink contours) crossed the penumbral 
boundary between t = 93.8 and t = 97.7 min and its proper 
motion outside the spot was -1.9 km s"^ Anomalous Stokes 
V profiles appear where one would locate the EC at ^ = 109.6 
min. We cannot calculate Doppler velocities from the Stokes V 
profiles observed in those pixels. However, the EC is easily de- 
tected in the velocity maps computed from the Doppler shifts of 
Stokes / (not shown). The evolution of EC E can be tracked until 
t = 109.6 min. Thus, EC E disappeared as a velocity structure 
1 1.9 min after crossing the spot boundary. 

The evolution of ECs M and N outside the spot is shown 
in Eigs. IB. II and IB. 21 EC M is only visible in two maps before 
it goes out of the EOV. EC N disappeared 2! I from the spot 
boundary. 
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Fig. 14. Same as Fig. [13] for the quantities inferred from Fe i 
630.25 nm. 
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Fig. 15. Maps of Stokes V zero-crossing velocity of Fei 
1565.28 nm as a function of time. Pink contours mark EC A. 
Orange contours enclose pixels with V630.25 and V1565.28 lower 
than -2.3 km s"^ and -2.7 km s"^ , respectively. The vertical 
lines mark a radial cut passing through the EC. Each tickmark in 
the y-Sixis represents V\ 



The phenomenological properties of ECs D, F, M, and N 
in the moat are summarized in Table [6l The behavior of type 
II ECs outside the spot is very similar: they quickly dissolve 
(-13.8 min) without reaching large radial distances from the 
outer penumbral boundary (~1V9). Inside the spot ECs show up 
as elongated structures, but type II ECs adopt a roundish shape 
when they enter the moat. Moreover, their proper motions are 
reduced relative to those they exhibit in the penumbra (compare 
Tables[6]and 1). 

As discussed by " Cabrera Solana et aL 1 (I2OO6I), ECs D and 
F became moving magnetic features (MMFs) in the sunspot 
moat. This suggests that type II ECs are the precursors 



Table 6. Phenomenological properties of type II ECs in the 
moat. The time and distance from the outer penumbral boundary 
where the ECs vanish as velocity structures are given in the sec- 
ond and third columns. The fourth column displays the time each 
of them survived in the quiet photosphere. The proper motions 
of the ECs in the moat are given in the fifth column. The last two 
columns display their maximum lengths and widths, corrected 
for LOS eff'ects, outside the penumbra. The parentheses indicate 
that EC M leaves the FOV before disappearing. 



EC 


^dis 


ddis 




^prop 


Length 


Width 




[min] 


n 


[min] 


[km s-i] 


[km] 


[km] 


D 


74.3 


1.8 


15.6 


0.6 + 0.3 


1300 


700 


F 


109.6 


1.8 


11.9 


1.9 + 0.4 


1000 


1200 


(M) 








3.5 + 0.0 


900 


400 


N 


57.8 


2.1 


14.0 


1.3 + 0.1 


900 


600 


Mean 




1.9 


13.8 


1.8 


1000 


700 



of at least some MMFs. The relation between MMFs and 
the magnetic field of the p enumb r a has been studie d by 
Sainz Dalda & Martinez PilleJ (l2QQ5h . iKubo et al.l (l2QQ7h . and 
iCabrera Solana et aL (.2008i) . 



13. Summary and conclusions 

We have presented a detailed description of the temporal evo- 
lution of the Evershed flow including, for the first time, its po- 
larization signatures. Our main findings can be summarized as 
follows: 

1. We identify fifteen Evershed clouds (ECs) that appear 
in the mid penumbra. They propagate with speeds of 
~2.6 km s"^ along penumbral filaments characterized by 
larger values of the linear polarization, linear-to-circular po- 
larization ratios, and Doppler velocities than their surround- 
ings. These filaments ca n be associated wi th the penumbral 
intra- spines described by lLites et alJ (1 19931) . 

2. ECs can be classified in two types: (a) ECs that disappear 
in the outer penumbra (type I); and (b) ECs that cross the 
visible border of the spot and enter the sunspot moat (type 
II). Most of the observed ECs belong to type 1. 

3. The motions of ECs in the penumbra are mainly radial. 

4. ECs have greater lengths (-1700 km) than widths 
(-700 km), and their sizes increase as they migrate outward. 

5. Our observations suggest the existence of diff'erent time 
scales associated with the appearance of the ECs. We find 
periodic behaviors of the order of 15 min and 40 min in two 
diff'erent intra- spines. 

6. ECs display larger Doppler velocities and L/V ratios in the 
infrared lines of the set. This implies that: (a) ECs are more 
conspicuous in deeper atmospheric layers; and (b) the use of 
infrared lines would be advantageous to detect them in vector 
magnetograms. 

7. ECs possess greater Doppler velocities, L/V ratios, and area 
asymmetries than the intra- spines along which they migrate 
outward. This might imply that the magnetic and dynamic 
configuration of the intra-spines are modified by the passage 
of ECs. Assumi ng that the velocity and m agnetic field vec- 
tors are aligned (iBellot Rubio et al .1120031) . the enhancement 
of Doppler velocity displayed by the ECs can only be pro- 
duced by an increase in the modulus of the velocity vector. 

8. The Doppler velocities observed in ECs are larger when the 
sunspot is farther from disk center, while their L/V ratios are 
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smaller. These results are consistent with the idea that the 
flow and the magnetic field are almost horizontal to the solar 
surface. 

9. ECs are brighter than intra- spines in the inner penumbra. 
The brightness excess is no longer seen in the mid and outer 
penumbra. 

10. The Doppler velocities of the ECs are systematically greater 
than the LOS projection of their proper motions, assumed 
horizontal to the solar surface. 

1 1 . The blueshifts, linear-to-circular polarization ratios, and area 
asymmetries of the ECs increase with radial distance. In the 
intra- spines, L/V and SA display the same behavior, while 
the velocity decreases from the middle to the outer penum- 
bra. 

12. The disappearance time scales of type I ECs are shorter than 
the cadence of ~ 4 min of the observations. 

1 3 . The outer penumbral boundary is modified by the passage of 
the ECs that reach the quiet photosphere. When a type II EC 
penetrates into the moat, the outer penumbral boundary does 
the same and the intensity decreases in the quiet photosphere. 

14. Four ECs are seen to cross the visible border of the spot. 
Once they enter the moat: (a) their propagation velocities are 
smaller than inside the penumbra (0.9 vs 2.6 km s"^); (b) 
they adopt a roundish shape; and (c) their blueshifts quickly 
decrease until they vanish. Type II ECs can be detected as 
velocity structures in the quiet photosphere for only about 
14 min after crossing the outer penumbral boundary. They 
disappear relatively close to the spot (-2''). 

15. Some type II ECs become mov ing magnetic features in the 
sunspot moat (Cabrer a Solana et al.n2006l) . but they do not 
show significant Doppler shifts. 

The line parameters derived from the polarization profiles 
give us valuable information about the phenomenological prop- 
erties of the ECs. However, the interpretation of these observa- 
tional parameters in terms of physical quantities is not straight- 
forward. For example, a larger value of the Doppler shift can be 
associated with a stronger Evershed flow, but also with a mag- 
netic field which is less inclined to the LOS. 

To remove these ambiguities, we have perform detailed 
Stokes inversions of the observed spectra. The results of the in- 
versions, presented in Paper II of this series, will allow us to 
investigate the nature of the EC phenomenon. 
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Appendix A: Removal of scattered light 

In Fig. lA.ll we present fits to the disk-center profiles of Fe i 
630.15 nm and 630.25 nm recorded on July 1, 2005 using Eq.[T] 
As can be seen, the two visible lines are well reproduced. The 
merit function of the fit (right panel) shows a minimum for 
K ^ 15% and cr ^ 6 pm. The same results are obtained for 
the profiles measured on June 30, 2005. 

Figure IA.2I shows an example of the eff'ect of the scattered 
light correction on the polarization spectra: the original Stokes U 
profile displays signatures of the telluric lines, but they disappear 
after the correction. 

Appendix B: Line parameter maps from the July 1 
observations 
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Fig. A.l. Comparison of the spectra at disk center of Fe i 630.15 nm (left panel) and Fe i 630.25 nm (center panel) as observed with 
FTS (thin solid) and POLIS (thick solid) for the 1 July dataset. The dashed lines are the best fits obtained by combining the FTS 
profiles with the fraction of scattered light and the resolving power of POLIS as explained in the text. The merit function of the fit 
is shown in the right panel. Note that the O2 telluric line is not considered in the fit. 




Fig. A.2. Stokes / (top) and Stokes U (bottom) profiles of the lines at 630 nm before (left) and after (right) correction for scattered 
light. The vertical lines indicate the positions of the O2 telluric lines. 
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Fig.B.l. Time evolution of line parameters in the center-side penumbra of AR 10781 on July 1, 2005. From top to bottom: 
Continuum intensity at 1565 nm, Stokes V zero-crossing velocity of Fei 1564.85 nm and Fei 630.25 nm, logarithm of the total 
linear polarization of Fei 1564.85 nm and Fei 630.25 nm. Color contours outline the ECs. The letters at the bottom of the third 
panel label each EC. Red lines indicate the inner and outer penumbral boundaries. Each tickmark in the j-axis represent V . The 
arrows marks the direction to disk center. 
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